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Abstract 

We present the study of exact inhomogeneous cosmological solu- 
tions to a four-dimensional low energy limit of string theory containing 
non-minimal interacting electromagnetic, dilaton and axion fields. We 
analyze Einstein-Rosen solutions of Einstein-Maxwell-dilaton-axion equa- 
tions and show, by explicitly taken the asymptotic limits, that they 
have asymptotically velocity-term dominated (AVTD) singularities. 

1 Introduction 

In the context of the so called string cosmologies there are several aspects of 
interest. One of them is the way they approach to the singularity as well as 
the asymptotic behaviour of spacetimes filled with coupled scalar and U(l) 
fields. Those matter fields arise naturally from low energy effective super- 
string theory [1] i.e. the Einstein-Maxwell-dilaton-axion (EMDA) system 
that contains the dilaton which is non-minimally coupled with the Maxwell 
and axion fields. The EMDA system has been discussed actively in the con- 
text of the black hole solutions and singularities, and possesses characteristic 
features different from the Einstein-Maxwell system [2] . 

In this paper we investigate the character of initial cosmological singular- 
ities in the presence of matter fields arising from the low energy superstring 
theory. 



Belinskii, Khalatnikov and Lifshitz (BKL) [3] conjectured that the dy- 
namics of nearby observers would decouple near the singularity for different 
spatial points. BKL also speculated that a generic singularity should exhibit 
such a local oscilatory behavior. There is a special case of the BKL conjec- 
ture called the asymptotically velocity-term- dominated (AVTD) singularity 
that is not described by an infinite sequence of Kasner epochs but by only 
one epoch. In this case no oscillatory approach to the singularity is observed. 
The characteristic feature of the AVTD singular behaviour is that all spatial 
derivative terms of the field equations are negligible sufficiently close to the 
singularity. 

Another aspect of interest is to investigate the influence of nonminimal ex- 
ponential coupling of the dilaton and axion to the [/(l)-field on the character 
of the singularities. As it has been shown in [3] and [4] the minimally cou- 
pled scalar field can suppress the oscilatory behaviour; on the other hand, it 
has been shown numerically that magnetic Gowdy spacetimes are not AVTD 
but mixmaster [5]. It is therefore important to study the behavior of models 
containing both, scalar and electromagnetic fields, in the limit t — > 0. 

Our aim is to verify BKL conjecture in non-vacuum and spatially inho- 
mogeneous cases. To this end we shall consider Einstein-Rosen spacetimes. 

The Einstein-Rosen models are important among other things because 
many of the non-perturbatively exact superstring backgrounds constructed 
from the gauged Wess-Zumino-Witten (WZW) models admit two Abelian 
isometries [6]. Hence, given our current knowledge of string cosmology and 
conformal field theory, the G2 Einstein-Rosen cosmologies derived within 
the context of the low energy effective action represent a set of models that 
is closely related to many exact string solutions in four dimensions. More- 
over, the Einstein- Rosen solutions include a number of spatially homogeneous 
Bianchi models as special cases and certain G2 models may be viewed as in- 
homogeneous generalizations of those Bianchi cosmologies. 

One subtle point in this treatment is that we are extrapolating asymp- 
totic behavior of solutions that are valid in the weak coupling regime to 
regions near the big-bang singularity, for instance. In the strong coupling 
regime one expects that higher-order corrections to the perturbative theory, 
as well as non-perturbative string effects, should become increasingly impor- 
tant. Hence, in this regime the qualitative behavior of these solutions may 
deviate somewhat from that of solutions derived from the full M or string 
theory. However, there are reasons to believe that G2 solutions should never- 
theless provide a generic description of cosmological models in the vicinity of 
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a singularity. A major incentive for this comes from the long standing con- 
jecture of Belinskii and Khalatnikov [3]. This states that on the approach 
to the cosmological singularity, the generalized Einstein-Rosen metrics may 
play the role of the leading-order approximation to the general solution of 
conventional Einstein gravity. 

This paper is structured as follows: We present the EMDA field equations 
and solutions in Sees. 2 and 3, respectively. In Sec. 3 we identify the previ- 
ous solution as an Einstein-Rosen spacetime and determine its kinematical 
parameters. We analyze the Raychaudhuri equation to verify the existence 
of singularity in Sec. 4. In Sec. 5 we perform the asymptotic analysis of the 
metric and the fields, and show that the approach to the singularity is of the 
form AVTD. In Sec. 6 we present the plane wave and cylindrically symmet- 
ric cases and show their approaching to t — 0. Finally, some conclusions are 
given in Sec. 7. 

2 EMDA Einstein-Rosen solution 

As we pointed out above, we shall consider solutions to the low energy ef- 
fective theory for the heterotic string that arises as a dimensional reduction 
and truncations of the string theory in four dimensions under the following 
considerations: compactifications of six of the ten dimensions and omission 
of the arising massless fields in the obtained heterotic structure; only U(l) 
charges are permitted. The dynamical equations of the resulting theory can 
be deduced from the action 

S = j d*x^{R - 2{d^f - - e^F 2 - «J^J>}, (1) 

where R is the scalar Riemann curvature, is the four-dimensional met- 
ric tensor, F^ v is the electromagnetic antisymmetric tensor field, F^ v its dual 
(F^ = —\\/—g€[ 1 vapF al3 ) 1 (p is the dilaton scalar field and k is the axion 
field dual to the three index antisymmetric tensor field H = — |e 4< ^ * dn. 
The Einstein frame is related to the string frame by a conformal transfor- 
mation: G^ v (s) = e^g^E), where s and E stand for string and Einstein, 
respectively. 

Previous work on these lines includes solution generating techniques which 
permit to construct exact inhomogeneous solutions to the equations of low 
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energy string theory containing nonminimally coupled dilaton and electro- 
magnetic fields (EMD) presented in [7]. Solutions to Einstein-dilaton-axion 
theory have been derived by applying the global symmetries of the string 
effective action to a generalized Einstein- Rosen metrics in [8]. An algorithm 
for generating families of inhomogeneous spacetimes with a massless scalar 
field was presented in [9]. 

In the Einstein-Maxwell-dilaton-axion (EMDA) theory, the dilaton 0, ax- 
ion k, gravitational g^ and electromagnetic Maxwell fields ought to fulfil 
the EMDA field equations derived from flTfl: 



In the next section we present EMDA solutions for metrics that possess 
two spacelike Killing vectors. The fundamental structural functions of these 
solutions are rational functions expressible as the ratio of polynomials of 
at most second degree in the coordinate variables. These solutions were 
obtained by straightforward integration of Eqs. (2)-(5) in [10]. 

2.1 Metric functions, fields and Weyl scalars 

The metric, endowed with two Killing vectors: d x and d y , is proposed in the 
form: 




(2) 



(3) 



(4) 
(5) 



ds 2 = ^dp 2 - ^-dq 2 + ^-(dx + Ndyf + Q(dx + Mdy) 2 , (6) 



The solutions are given by 



4 



P = ep 2 + 2np + a, 

Q = eq 2 + 2mq — a, 

M = up 2 + 2bp, 

N = -vq 2 + 2(3q, (7) 

together with e = 1, 0, — 1, v = 1, 0, —1, A = M — N . Being m, n, b, (3, a 
constants. There are distinct families of solutions depending on the value of 
the kinematical constant e = 1,0,-1. The corresponding matter fields are 
given by 

2<b, \ uj(p 2 + q 2 ) 
e*{P,q) = ^ , (8) 

2(bq + f3p) 
w{p + q ) 

a.m = do) 

where k and u are axion and dilaton parameters, respectively; while is 
the electromagnetic potential and go and go are constants related to the elec- 
tromagnetic field. To fulfill the EMDA equations, the constants m, n, b, (3, a, g , q , k , u, e, v 
are restricted to satisfy the following algebraic conditions: 



v 2 ql = 2uP(Pe + mis), g (3 - q b = 0, (12) 
v 2 gl = 2ub{be-nu), nf3 + mb = 0, (13) 

This constrictions arise in the process of solving the field equations with 
the metric functions proposed as polynomials. As a result, not all the con- 
stants can be considered as free parameters; from the eleven constants only 
five of them are free parameters. 

The nonvanishing Weyl scalars are given by; 



= ^ = -^r(& 2 + A (14) 
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6A 3 ip 2 = — {Qv{nv — eb)[p(p 2 — 3q 2 ) + iq(q 2 — 3p 2 )} 
—6u(mu + e(3)[q{q 2 — 3p 2 ) — ip(p 2 — 3q 2 )] 
+2[2e(6 2 + (3 2 ) + 3v(m/3 - nb)} (q 2 - p 2 ) 
+12ipq[(be — nv)b + {vm + e/3)(3] 
+A(nb - m(3)(bp + Pq) + 4a(/3 2 + 6 2 )}, (15) 

Since the three Weyl scalars ipi,ip2, ip3 are nonzero, the solution is of type 
G in the Petrov classification. From the above expressions it can be seen that 
it is the presence of the constants b and j3 (non-diagonal terms in (JHJ)) that 
gives such structure, otherwise being a type D metric without axionic nor 
dilatonic fields; these fields become constant when b = = j3. The invariants 
that can be constructed from the Weyl scalars are = 6(-^2) 2 — 8ipiip3 
and = 12^1^3 - 6(^ 2 ) 3 - It can also be noted that such invariants 
will diverge provided the Weyl scalars do. This is the case if the function 
A = M — iV is zero, then essential singularities arise at those points. 



3 Einstein-Rosen spacetimes 

The line element (JHJ) can be transformed into the generalized Einstein-Rosen 
form: 

ds 2 = e f (dz 2 - dt 2 ) + lah dx a dx h . (16) 

where x a = x,y. All components are independent of the spatial coordi- 
nates (x,y); f = f(z,t) determines the longitudinal part of the gravita- 
tional field. The metric of the surfaces of transitivity is 7 a b and the gradient 
— d^{det^ ab ) 1 / 2 determines the local behaviour of the model. If is 
globally spacelike or null, the solutions represent cylindrical and plane grav- 
itational waves, respectively. The cosmological models are characterized by 
a timelike or when the sign of K^K^ changes. 

In terms of the metric functions A, P and Q of (JHJ) the Einstein- Rosen 
metric is 

ds 2 = A(dz 2 - dt 2 ) + ^{ X (dx + Ndyf + x~\dx + Mdy) 2 } (17) 

where \ = \JPjQ and G = y/PQ. The appropriate coordinate transfor- 
mation (p, q) i — ► (z, t) that leads from (JHJ) to (JT7J) depends on the values of 
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the constants e,n,m,a. In Tabled some cases are shown. 



e 


P 


q 


conditions 


1 


\Ja — n 2 sinh z — n 


\Ja + m 2 cosh t — m 


a — n 2 > a + m 2 > 


1 


e z /2 - n 


e t /2 — m 


a = n 2 = —m 2 = 


-1 


\J a. + n 2 sin z + n 


\/ m 2 — a sin t + m 


a + n 2 > m 2 — a > 


1 


e z /2 - n 


\J a + m 2 cosh t — m 


a = n 2 


1 


\^n 2 + m 2 cosh z — n 


e t /2 — m 


a = — m 2 





z' 2 n 2 —a 


t 2 +a 




2n 


2m 





Table 1: Coordinate transformations (p, q) — > to pass from the line element 
© to the Einstein- Rosen form (|17|) . 

The calculus of the gradient of the transitivity surface shows that the 
three possibilities of interpreting the solution (fTTj) are available as particular 
cases. We are interested in particular in the cosmological model obtained 
with the transformation (p, q) i— > (z, t) given by 

p *—>■ V a — n 2 sinh z — n, 

q i— > V a + m 2 cosh t — m. (18) 
This transformation in the line element (|17|) gives 

Ads 2 = A 2 (dz 2 - rft 2 ) + (a - n 2 ) cosh 2 z(dx + Ndyf 

+ (a + m 2 ) sinh 2 t(dx + Mdy) 2 , (19) 

where A = M(z) — N(t) is a monotonically increasing function on z and t, 
given by 



A = v{y a — n 2 sinh z — n) 2 + 26(V a — n 2 sinh z — n) 

+i/(v a + m 2 cosht — m) 2 — 2/3(v a + m 2 cosht — m). (20) 

In our case A ^ over the whole spacetime and it implies that the invari- 
ants are regular as well. The fact that the invariants are regular and that the 
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spacetime has a nonvanishing acceleration, as it is shown in the next section, 
might be an indicative of absence of singularity. A nonvanishing acceler- 
ation may create a gradient of pressure which acts opposing gravitational 
attraction and this fact allows to avoid the focusing of the congruence in the 
Raychaudhuri equation. A singularity-free spacetime due to its acceleration 
was given in [11]. In the next subsection we characterize the cosmological 
model using kinematical parameters. 



3.1 Kinematical parameters 

The kinematical characteristics of the Einstein-Rosen spacetime (j!7|) are de- 
termined with respect to the timelike vector u a = 7^^", u a u a = — 1. The 
acceleration, expansion, deceleration and shear are given, respectively, by: 



ill 







M 
2A 2 ' 

QA-NQ 
2A 3 / 2 Q ' 

q = -. —. {8N 2 Q 2 + 5A 2 Q 2 

(QA - NQ) 2 

-AQNQ + QAQ{QN - AQ)}, 
,2 Q 2 4QiV iV 2 (3P + 4Q) 

6cT = a^ + a^q + — a^q — ' (21) 

being null the rotation or vorticity, uj^ = 0. For the cosmological model the 
kinematical parameters behave asymptotically as follows 



e(t -»• 0) -> 00, e(t -> 00) -> 0. (22) 

a(t -4 0) -> 00, tr(t 00) ->• 0. (23) 
g(t — > 0) — > const, q(t — > 00) — > const. (24) 

This behaviour is illustrated in Fig^ The shear evolves towards isotropy 
for large times, a 2 —>■ 0. The expansion, being infinite at the origin of time, it 
decreases as time passes. The deceleration is positive all the time and tends 
to be constant, this indicates that the spacetime expands in non accelerated 
way. Both, shear and expansion diverge at t — 0, indicating the presence of 
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Figure 1: Illustration of the behaviour of kinematical parameters for the cos- 
mological solution. Deceleration parameter, q, shear, a, and the expansion, 
O, are shwown. For this plot the constants have been fixed to m = n = 2, 
b = —(3 = 1, a = 5 and v = 1. 

a singularity; the analysis of the Raychaudhuri equation in the next section 
confirms that the spacetime has a spacelike singularity at t — 0. 

4 Raychaudhuri equation 

According to Raychaudhuri, the equation that governs the rate of change of 
expansion of timelike congruences, O, is 



where u a is a tangent vector and A is the affine parameter of the congruence, 
Rob is the Ricci tensor, u is the rotation or vorticity, o is the shear and ii a is 
the acceleration of the congruence. 

From Eq. ()25|) it can be seen that the expansion of a timelike geodesic 
congruence with zero vorticity will monotonically decrease along a geodesic if, 
for any timelike vector V a , R a bV a V b > 0, i.e. if the strong energy condition 
(SEC) is hold. The term involving the Ricci tensor, R a b, in Eq. (|23j) induces 
contraction of the geodesic lines, indicating that the focusing of neighbour- 
ing geodesies is unavoidable if SEC is fulfilled since the other terms on the 
right hand side of (|25j) are also negative, with exception of the acceleration; 




:a 



(25) 
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however, as we shall see the positive term of the acceleration is defeated by 
the term R ab u a u b . Substituting from ( [2T|) in Eq. (|2*5J) we obtain: 



+ 4^? {M T ~ 3Ar § + 2M " + 2jV) ' (26) 

Moreover, from the expression for B = dQ/dX, Eq. (|25|). it can be seen 
that in fact it diverges as t approaches the origin provided the metric function 
Q be divergent, 

2A X Q Q 2 2QA 2A 2 A 1 ; 

In the cosmological case, Q = (a + m 2 )sinh 2 t, therefore as t — > 0, 
8 « -1/ sinh 2 t — > — oo, revealing the existence of a singularity as t — > 0. 
This behaviour can be confirmed by writting the derivatives of the coordi- 
nates with respect to an affine parameter and checking that their coefficients 
diverge. By showing that the first derivatives are unbounded one concludes 
that the corresponding geodesic curves t(r), x(r), y(r), z(r), are uncomplete 
at t = [12]. 



5 Asymptotic Behaviour 



In the limit t — > the fields become functions that depend on z but otherwise 
are constant, 



e 2 *(t -> 0) 



u(A 2 + (Va + m 2 — m) 
A 



. n , 2b{ Va + m 2 -m)+ 2(3 A 

K{t — ► 0) — > , h «o, 

u;(v4 2 + (Va + wi 2 — m) 2 ) 



A x (t -> 0) 



g> ^4 - g (VaTtf? - m) 
A 



d(Vtt + ff?-m)M + So(i/aTn? - m)) 
A y (t — > OJ — >■ — , 
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In relation to the asymptotic behaviour as t — > oo, both scalar fields, 
axion and dilaton tend to constant values as t approaches infinity, therefore 
the fields decouple for large times, 





2 v 



n{t — > oo) — > k , 
A x (t ^ oo) -> 0, 
Ay(t oo) -> v4# , 



(29) 



There exists a lower bound on the value of the dilaton field and this implies 
the existence of a lower, non- vanishing bound on the string coupling which, in 
the context of M theory, in turn implies the existence of a lower bound on the 
radius of the eleventh dimension [13]. Moreover, in the limits that the volume 
of the transverse space measured in the string frame becomes vanishingly 
small or arbitrarily large the axion field is constant, since F =detr a f, = e^PQ, 



lim T = — (a — n 2 )(a + m 2 )(cosh 2 z)e 2t — > oo, k — > k . (30) 

t— +CJO 4lU 

Thus the two-form potential effectively decouples from the field equations 
in these limits. 

As was mentioned at the begining, one of the most interesting aspects of 
analyzing cosmological models in to dilucidate how the spacetime behaves 
near the initial singularity and to determine if this is an oscilatory approach 
or a AVTD behaviour. 

Narita, Torii and Maeda [14] studied the influence of the exponential 
dilaton-electromagnetic coupling on the character of initial singularities. They 
considered EMDA system with a T 3 Gowdy cosmology and using the Fushian 
algorithm they showed that those spacetimes have in general asymptotic 
velocity-term dominated singularities. Their results mean that the exponen- 
tial coupling of the dilaton to the Maxwell field does not change the nature 
of the singularity. In what follows we arrive to the same conclusion for the 
Einstein- Rosen spaces by taking directly the limit t — > in the metric ([19|. 

Taking the limit t —>■ with z =constant, in the metric functions of (|T9~|) 
for the studied case we obtain the line element 



limT 



(a - n 2 ) (a + m 2 ) (cosh 2 z)t 2 



k — > const, 



11 



P(t\ (rv 4- m 2 )t 2 

ds 2 = k(dz 2 - dt 2 ) + -^-(dx + Ndy) 2 + 1 , ' (dx + Mdy) 2 , (31) 

k k 

where k =const. and the metric functions are ( remind that z has been fixed 
to constant) 



P( z ) = ( a -n 2 )cosh 2 (z), 

M(z) = v\\J a — n 2 sinh(,2) — n} 2 + 2b[V a — n 2 sinh(,2) — n], 

N = -v[Va + m 2 - m] 2 + 2/3{Va + m 2 - m). (32) 

The line element ()3ip corresponds to a Kasner metric with pi = p 3 = 
and pi = 1. This establishes the AVTD behaviour of this model, extending 
the conclusion by Narita et al to Einstein- Rosen spaces. 



6 Plane-symmetric Wave and Cylindrical space- 
time 

To establish a comparison of the behavior near the singularity, we take the 
limit as t — > in the plane-symmetric wave spacetime and the cylindrical 
one associated to the metric (fTTjl . As quoted in Sec. 3 the local behavior of 
the Einstein- Rosen spacetime has to do with the character of the gradient of 
the transitivity surface. The spacetime (|T7|) can be interpreted as a plane- 
symmetric wave space by using the coordinate transformation 

p i— > exp z/2, q^expt/2, (33) 
then the metric (JT7jl takes the form 

2z 2t 

ds 2 = A(dz 2 - dt 2 ) + —(dx + Ndy) 2 + —(dx + Mdy) 2 . (34) 

This metric behaves as a conformally plane spacetime when t — » 0: 

A 2 1 
ds 2 = k(dz 2 - dt 2 ) + —(dx + Nxdy) 2 + — (dx + MM 2 . (35) 

4k 4k 

where A = exp 2z, Mi = v 2 exp2z/4+f3 exp z are constants for z = const, 
and Ni = —u 2 /4 + b. The metric is completely regular at t — > 0; besides, 
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neither 6 nor a diverge at t — 0, in contrast with the cosmological case. 
Investigating the tendency of the change in expansion along a congruence, 
O, we found that it is finite for all the time range. However the possibility 
exists of a fine tunning of the axion and dilaton parameters, b and j3, that 
make A = and in that case the spacetime becomes singular. 

The cylindrically symmetric case arises if instead we make the coordinate 
transformation, 

p i— > \J m 2 + n 2 cosh z — n, q h- > expt/2 — m. (36) 

In this case the norm of the gradient of the transitivity surface is always 
positive, G a G a = e 2t (m 2 + n 2 )/(AA) and the spacetime can be interpreted as 
a cylindrical symmetric one; the form of the metric is: 

ds 2 = A(dz 2 - dt 2 ) + K + ^) sinh2z (rfx + Ndy) 2 + ^(dx + Mdyf. (37) 
In this case the singularity at z = is timelike. 

7 Conclusions 

We studied AVTD cosmological Einstein-Rosen spacetimes with EMDA fields. 
The matter fields show a tendency to decoupling for large times. Near the 
singularity the dynamics at different spatial points decouples and the metric 
has a spatially varying Kasner form. We conclude that the nonminimally cou- 
pling of the dilaton and axion to the Maxwell field does not change the nature 
of the singularity. In this sense we have extended to generalized Einstein- 
Rosen spaces the result, founded in Gowdy spaces by NTM, that exponential 
coupling of the scalar field does not necessarily lead to Mixmaster behavior. 

Further analysis on solutions with plane-wave and cylindrical symmetry 
might be of interest as roughly is shown in Sec. 6. The asymptotics of the 
fields resemble the cosmological case, however, each model approaches t — 
in a very different way. 
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